Arrested States formed on Quenching Spin Chains with Competing Interactions and 

Conserved Dynamics 
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We study the effects of rapidfy coofing to T = a spin chain with conserved dynamics and 
competing interactions. Depending on the degree of competition, the system is found to get arrested 
in different kinds of metastable states. The most interesting of these has an inhomogeneous mixture 
of interspersed active and quiescent regions. In this state, the steady-state autocorrelation function 
decays as a stretched exponential ~ exp(— (£/ t °) 3 )j an< i there is a two-step relaxation to equilibrium 
when the temperature is raised slightly. 



PACS numbers: 05.70.Ln, 05.40. +j, 81.40.Ef 

When a system at a high temperature is cooled rapidly 
to low temperatures, it may not be able to reach an equi- 
librium state in an experimentally realizable time. In- 
stead it may reach a long-lived arrested state, often with 
a degree of spatial disorder, which ultimately relaxes to- 
wards equilibrium over very much longer time scales. An 
intriguing possibility is that the disorder induced by the 
kinetics may be strong enough that the system has widely 
different levels of dynamical activity in distinct regions 
of space. Dynamically heterogeneous states are found to 
arise, for instance, in a glass- forming liquid [|l|. From the 
theoretical point of view it is important to ask: Are there 
simple models in which dynamically heterogeneous states 
arise naturally? Can their formation and properties be 
understood in microscopic terms? Finally how do such 
states decay, and how is equilibrium approached? 

We address these questions by studying nonequilib- 
rium quenches to T = 0, in simple lattice models. Fol- 
lowing such quenches, the system may get arrested in a 
metastable state instead of reaching the ground state. A 
useful way to characterize the resulting arrested state is 
to ask whether or not there is any dynamical activity in 
it. A quiescent arrested state is one in which the system 
settles into a single configuration, and degrees of free- 
dom are frozen. Another possibility is that the arrested 
state may involve a large number of configurations which 
are dynamically accessible from each other; in that case 
the system is dynamically active. Interestingly, all these 
possibilities are realized in quenches of a family of sim- 
ple models, namely Ising chains with different degrees 
of competition and conservation laws. In the absence 
of competition, the system approaches the ground state 
if the dynamics is nonconserving |2|], while it reaches a 
uiescent arrested state under spin-conserving dynamics 
On the other hand, a system with competing interac- 
tions, evolving under nonconserved dynamics, has been 
shown to exhibit an active arrested state Q. This natu- 
rally leads to the question: Are new features brought in if 
both conservation and competition are present? In this 
Letter, we study the effects of quenching the simplest 
model which incorporates both these features, namely 



an Ising model with competing first and second neigh- 
bour interactions, evolving through a dynamics with a 
single conservation law. Despite its simplicity, the model 
shows interesting transitions in the character of the ar- 
rested states as the degree of competition is varied (Fig. 
1). The most interesting of these arrested states, reached 
for strong enough competition, is of a qualitatively new 
type: it has active and quiescent regions interspersed in a 
disordered fashion throughout the system. This inhomo- 
geneous quiescent and active (IQA) state, has nontrivial 
dynamical properties. In the T = steady state the 
autocorrelation function has a stretched exponential de- 
cay. Further if the temperature is raised slightly, there is 
a two-step relaxation: the IQA state relaxes to equilib- 
rium via an intermediate long-lived intermediate energy 
state. We are able to quantitatively understand many 
of these unusual features, often found in glassy systems, 
within this simple model. 

The equilibrium phases and transitions of the Axial 
Next Nearest Neighbour Ising (ANNNI) model have been 
well studied and characterized [pi. However, its nonequi- 
librium properties remain relatively unexplored even in 1- 
d, except for a few studies. The first such study explored 
arrested states obtained by quenching across T — phase 
boundaries of an extended ANNNI model appropriate to 
polytypes Q. More recently, time-dependent coarsening 
induced by nonconserved dynamics has also been stud- 
ied by quenching the system across phase boundaries at 
T = @], and also from T = oo to T = [g. Here we 
explore the interplay of competing interactions with con- 
servation laws in the dynamics. Our principal result is 
the identification and characterization of an IQA arrested 
state in this system. 

Consider an Ising chain with spin variables {si} de- 
scribed by the Hamiltonian 



h = -Ji)^sjSi + i + ^y^s 



i+2- 



(1) 



An antiferromagnetic next neighbour coupling (J2 > 0) 
competes with the nearest neighbour coupling J\ which 
may be of either sign. In what follows, we shall assume 
J 1 > (ferromagnetic coupling) and define ji — J%/ J\ 



as a measure of the strength of competition. The equilib- 
rium ground state shows a transition from the ferromag- 
netic state tilt ■•• f° r 3i < 0-5, to the antiphase state 
tfll ••■ for j2 > 0.5. The point j'2 =0.5 is a multiphase 
point, at which the number of ground states (all configu- 
rations with no single spins) is exponentially large in sys- 
tem size L m . We are interested in the effect of a quench 
from an infinite-temperature random configuration of the 
system to T — 0. We use a double-spin-flip dynamics 
(DSFD), in which an adjacent pair of randomly chosen 
parallel spins is flipped (... tt •■• — > •■■ II •■■)■ Flips 
are attempted at unit rate, and are allowed only if the 
energy is not raised (AU<0). Evidently, the DSFD con- 
serves the difference M = M\ — M.% of the two sublattice 
magnetizations M\ and M^ . The dynamics thus involves 
a single conservation law. The DSFD maps onto the 
well-known Kawasaki spin exchange dynamics through a 
sublattice mapping, in which every spin on one of the 
two sublattices is inverted and the sign of the nearest 
neighbour coupling J\ is reversed. In the mapped model 
with Kawasaki dynamics, M is the total conserved mag- 
netization. We will use the DSFD, rather than Kawasaki 
description; results can be translated readily. The DSFD 
can be looked upon as an extension of the single-spin-flip 
Glauber dynamics to multiple-spin flips at a time. Multi- 
spin moves arise in physical contexts such as stacking dy- 
namics in the 3C-6H transition in SiC M and deposition- 
evaporation dynamics of bunches of particles || . 

The energy non-raising condition, which is a conse- 
quence of a T = quench, imposes local constraints on 
whether or not a pair of chosen spins can actually be 
flipped; these constraints are a function of ji- The nor- 
malized energy changes Ae=AE/ Ji involved in flipping 
a pair ff to || depend on the environments of the pair, 
and are given in Eq|2 below. There are six distinct lo- 
cal environments; the other unlisted environments are 
related to these by reflection symmetries. 
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The reverse of move (a) in (H) will be referred to as (a) , 
and similarly for the others. 

Evidently the dynamics is identical for all values of j'2 
for which the same set of moves are allowed. Moves (b), 
(e), (e), (/) and (/) are allowed (i.e. Ae < 0) for all 
]2- As ]2 is varied, Ae changes sign for moves (a), (c) 
and (d) . Each such change causes a change in the nature 
of the arrested state. While (c) is an allowed move for 
22 < 0, (c) becomes an allowed move for y\ > 0. Sim- 
ilarly, across ji — 0.5 and ji = 1, the allowed moves 
change from (d) to (J) and (a) to (a) , respectively. Thus 
there are distinct regions of dynamical activity along the 



j 2 axis: (i) j 2 e(-oo,0), (ii) j 2 e[0,0.5), (iii) j 2 = 0.5, (iv) 
j2£(0.5, 1], and (v) j2G(l,oo); see Figure 1. In region 
(v), i.e. for strong competition, the system reaches an 
IQA arrested state. 

We used Monte Carlo simulation to study the arrested 
steady states that are reached under Ae < DSFD start- 
ing from a random initial configuration corresponding to 
T = oo. We studied the approach to the arrested states, 
the dynamical behavior in these states, and finally the 
relaxation from these states to equilibrium at low but fi- 
nite temperatures. The approach to, and the decay from, 
the steady state was monitored by following the decay of 
the energy in time. Further, in cases (iii) and (v), we 
studied the dynamical behavior of the steady state by 
monitoring the spin-spin autocorrelation function 



C W = ^X^.fo&Mto + *)> - («i(*o)) S 



(3) 



where t is the number of Monte Carlo steps per spin 
and (...) denotes an average over t . We allowed for an 
explicit dependence of averages on the space location i, 
as arrested states need not be translationally invariant. 
Only at the multiphase point (iii) is the ground state 
reached on quenching; in the other four regions of j-i dis- 
cussed above (Fig. 1), the steady states are arrested. 

Before discussing the IQA state in detail, we sketch 
some features of the states in the other four regions. 

(i) The arrested state is quiescent. It consists of ferro- 
magnetic patches separated by clusters of frozen domain 
walls, e.g. ... TTTTITITTT ■■■■ It is qualitatively similar to 
the arrested state obtained in [[|, with only first neigh- 
bour interactions. 

(ii) The steady state has a number of diffusing domain 
walls separating ferromagnetically aligned patches (Fig. 
1). Though it resembles the active arrested states found 
in H, there is an important difference. The level of ac- 
tivity is much lower in our case, as the number of walls 
increases as ~ L? as opposed to ~ L in [Q. 

(iii) At the multiphase point there is a large degree of 
activity (Fig. 1), because the Ae = moves (d), (f) 
and their reverses carry the system through a subspace 
of ground state configurations labeled by a given value 
of M, The autocorrelation function C(t) ~ £~ 5 at long 
times, as in the unconstrained DFSD M. 

(iv) The steady state has alternating opposite-spin clus- 
ters of 2 or 3 spins, e.g. ... UITTIIITTT ■••> an( i is quies- 
cent. A single cluster of 4 or 5 spins may remain in the 
steady state and diffuse through a quiescent background 
(Fig. 1). 

In region (v), an IQA arrested state with alternating 
quiescent and active stretches, is reached. A segment of 
a typical configuration is depicted below. 

■• [nun] ITIIITITITTTITI [MMMMMMtMM]- 

Quiescent Active Quiescent 



Each active region has parallel-spin triplets in a back- 
ground of alternating single spins, while the quiescent 
portions predominantly resemble the arrested state of re- 
gion (iv) . Crucial to the coexistence of active and quies- 
cent regions is the existence of stable walls at the bound- 
aries of quiescent regions. These consist of left bound- 
aries fl~.IJ.JJ. or JJ-frfr and right boundaries JJ-JJ-ff or frfrJJ-, and 
are stable as moves (a) and (c) are energy-raising for 
J2 > 1. The numbers of quiescent and active (q and a 
respectively) regions of size £ are found numerically to 
decay as exp{—\t) with \ q ~ 0.05 and A a — 0.25. 

We now turn to the dynamical properties of the IQA 

state. The autocorrelation junction in the steady state 

i 
decays as a stretched exponential ~ exp(— (t/r ) 3 ) (Fig. 

2). Interestingly, the dynamical behavior of the IQA 
state can be related to the well-known symmetric ex- 
clusion process (SEP) of particles on a line |l(j]. This 
can be understood as follows. Each spin triplet in an ac- 
tive stretch can move by one unit right or left, under the 
DSFD move (e) (... |TI TTT IT - «— » - IT III TIT ■■•)• 
There is a hard-core repulsion between triplets as move 
(c) is disallowed. The dynamics within a single active re- 
gion is thus precisely that of a symmetric exclusion pro- 
cess of hard triplets on a lattice, where the single spins 
can be viewed as holes. The autocorrelation function 
C^(p, t) averaged over spins, of an active stretch of length 
I with p£ triplets (where £ = £ + 2 and 1/t < p < 1/3), 
is thus governed by the diffusion of these hard triplets; 
triplets extend over an extra lattice unit at both bound- 
aries of an active stretch making its effective length £. 
Hence, we expect Cg{p, t) to decay as t~? for times t less 
than a cutoff time rt(p), and as exp(—t / Tt(p)) thereafter. 
Further, Ti(p) can be found by noting an exact mapping 
of every configuration of this problem to a correspond- 
ing configuration of the SEP. Under the mapping, every 
triplet is replaced by a single particle, while a single spin 
maps onto a hole (TITIIITITTTIIITI — > ooo»oo»»oo). 
The mapped chain has a reduced length £' = £(1 — 2p). 
The stochastic W-matrices for the two processes are 
same, as there is a 1-1 correspondence between config- 
urations and moves. This implies that the eigenvalue 
spectra of the W-matrices in the two problems are the 
same, and in particular, the gap Ai to the first excited 
state is the same. The inverse of the gap is just the cut- 
off time T£, and so the above equality implies Tg = r' e ,. 
For the exclusion process, with free boundary conditions 
r'f, = 2£' /it 2 for large £' (the diffusion constant = ^) 
0, and hence r e (p) = 2£ 2 {\- 2p) 2 /n 2 . 

The autocorrelation function Cjqa (t) of the IQA state 
can be expressed in terms of a sum over active stretches: 



C IQA {t) = Y J Pi{p)c i { P ,t) 



(4) 



i'.f> 



where Pi(p) is the probability of finding an active stretch 
of length £ and density p of triplets. Even without explic- 
itly determining Pji(p), we can derive bounds on CiQA{t) 



using Pi — S p P|(p) ~ exp(— \ a £). As p varies across its 
range, the cutoff time T£ varies between the two limits 
72(0) = 2£ 2 /tt 2 (for a single triplet) and r £ (±) = 2£ 2 /9tt 2 
(for a single hole hopping over 3 lattice units). For each 
of these limits n(p*), the sum in Eq.(0) is dominated 
at long times by the term with the saddle point value 

£* = {tTT 2 /X a (l - 2p*) 2 y. The bounds imply that C IQA 

i 
has a stretched exponential form ~ exp(— (t / t ) 3 ) , with 

8/2437T 2 < T X 2 a < 8/27tt 2 . The numerically determined 

values r ~ 0.08 and A a ~ 0.25 are consistent with these 

bounds (Fig. 2). 

The dynamics of approach to the IQA state, starting 
from a random initial configuration, is also interesting. 
From numerical simulations the energy is found to de- 
cay as ~ exp(—(t/t )3) (see inset in Fig. 2). This is 
associated primarily with the fall in the number N& of 4- 
spin clusters which diffuse through ground state stretches 
••■ TTII ■•• until they dissociate when they encounter 
'traps' in the form of a single spin or a triplet; e.g. 
... TTITTTT .- - ». TTIIITT -■ The typical time for 
a 4-cluster to diffuse over a length I before encounter- 
ing such a trap is l 2 /D, implying that N^l) decays as 
exp(—Dt/l 2 ). Further, the stretch lengths are distributed 
exponentially (~ exp(-Xl)), so that the average of N^(l) 
over I is dominated by a saddle point value I* = (2Dt/X) 5 
at large t. This argument is reminiscent of that in ]l2J] 
and implies a stretched exponential form for the decay. 

Finally, let us discuss the relaxation to equilibrium from 
the IQA arrested state. The IQA state has regions with 
two distinct types of excitations, namely active patches 
with mostly single spins (and occasional mobile triplets) 
and quiescent patches with mostly triplets (and occa- 
sional frozen single spins). If T is raised to a small value, 
the system relaxes to an equilibrium state close to the 
... || J, | ... ground state, by annealing out both the single- 
spin and triplet excitations. Figure 3 shows the subse- 
quent variation of energy with time. There is a relatively 
rapid approach to a second metastable state, evidenced 
by a long plateau (Fig. 3), followed by an eventual ap- 
proach to equilibrium. This can be understood as follows. 
For finite T, the energy raising moves (a), (b), (c) and 
(d) are allowed, with probabilities uik ~ exp(—A.ek/T), k 
— a,b,c,d; the associated time scales are ~ 1/oJk- Moves 
(a) and (c) are instrumental in annealing out the two 
different kinds of local excitations (isolated spins and 
triplets respectively). They act on time scales which are 
widely different (ui^ 1 /ui^ 1 ~ exp(4Ji/T)) leading to the 
plateau. For small t , i.e. t ~ w" 1 , only move (a) is 
effective, which destabilizes active-quiescent boundaries 
and creates 5-spin clusters. Single spins diffuse out of 
active stretches and annihilate on meeting 5-clusters e.g. 
IITIIIII— ►UTTTIU- After the single spins anneal out, 
the system reaches a metastable state with clusters of 
length 2 and 3, much like the arrested state in region 
(iv). This continues till t > w^ 1 , when triplets begin 
to decay. To leading order in low T, the predominant 
decay channel involves the following steps: (?) the con- 



version of a triplet to a single spin and a 4-spin cluster 
(at rate u> c ), (ii) the production of a single spin when the 
4-cluster meets the nearest triplet, and (Hi) the fast dif- 
fusion of single spins till they meet triplets or single spins 
at separation In (n = odd), whereupon they annihilate 
e.g. UTUTTT^Umm^Umm. Process (Hi) is a 
variant of the single species diffusion-annihilation process 
||13| , implying a power-law (~ t~z) decay for the energy. 

To summarise, a simple understanding of the dynamics 
of the IQA state can be achieved in terms of diffusing ex- 
citations; the nature of approach, steady-state autocorre- 
lation function and decay of the state involve variants of 
the diffusion problem. For instance, the approach to the 
IQA state involves diffusion in the presence of randomly 
placed traps, while the autocorrelation function involves 
the consequences of confinement of diffusing excitations 
in active stretches of random lengths. In both cases, an 
average over the dynamically generated randomness re- 
sults in a stretched exponential decay. The two distinct 
time scales for relaxation from the IQA state arise from 
the different activation rates for the two types of diffusing 
excitations. Diffusion-limited annihilation of the second 
type governs the power-law decay towards equilibrium. 

We conclude by pointing out that IQA arrested states 
occur in several other situations, for instance, with anti- 
ferromagnetic nearest- neighbour coupling (J\ < 0), and 
also under a quench from a quiescent arrested state in re- 
gion (i). Further an IQA state is found in quenches of an 
extended ANNNI model relevant to polytype transitions 
|| . This model is richer, and shows variability in the 
microscopic nature of activity in IQA arrested states for 
different parameter values JlJ] . Interestingly, despite this 
variation, the dynamical behavior remains of the same 



form — a general consequence of the diffusion-based de- 
scription given above. 
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Figure Captions 

Figure 1: Space-time depiction of the activity (shown white) in different arrested states. Different regimes: (i) 

Quiescent (ii) Mobile domain walls (iii) Active (iv) Essentially quiescent (v) IQA state. 

Figure 2: The autocorrelation function in the IQA state with L — 12000 and 10 6 histories. The dotted curves are 

the bounds discussed in the text. Inset: Decay of energy excess (in units of J±) over the IQA value, starting from a 

random state. 

Figure 3: Energy per site e(t) (•) in units of J%, measured from the ground state value, when the IQA state is taken 

to T ~ 0.29Ji. We used u> a = 1, w c = 0.9 x 10~ 6 , L = 1200 and 12 initial conditions. Also shown are the fraction of 

single spins (A) and triplets (□). 
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